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Abstract 

We consider the heavy quark-antiquark potential in the AdS / QCD correspondence, focusing in 
particular on a recently calculated AdS-like metric deformed by back-reaction effects. We find that 
tuning the long-distance behaviour of the potential leads to a discrepancy at small distances, and discuss 
how to better constrain AdS / QCD geometries. A systematic comparison of various geometries is 
presented, based on goodness of fit to lattice data in the quenched approximation. The back-reacted 
geometry is seen to be phenomenologically favoured over an alternative geometry with the same number 
of parameters, although it does not perform as well as some other geometries. 



1 Introduction 

Regardless of whether string theory provides a fundamental description of nature, there is mounting evi- 
dence in favour of its being a useful tool for understanding strongly coupled gauge theories. The AdS / CFT 
correspondence of Maldacena [I] conjectures a mathematical equivalence between the low energy limit of 
type IIB string theory on AdS 5 <g)S 5 and a conformal field theory, Af = 4 U(N) Super- Yang-Mills (SYM) for 
large N, on the boundary of this space. Since this initial conjecture, various generalisations of the duality 
principle have been considered involving other gauge theories and their proposed gravity duals. Crucially, a 
weakly coupled string theory corresponds to a strongly coupled gauge theory, and it is this apparent fact that 
generates much phenomenological interest in such correspondences. The hope is that a geometry might be 
found whose dual field theory mimics Quantum Chromodynamics (QCD), thus facilitating the use of string 
theory to ascertain the properties of strong interactions. There are two main approaches. Firstly, one can 
develop consistent string theories and try to ascertain the properties of their field theory duals, with the hope 
of edging closer to QCD. Such theories are ten dimensional and must then be compactified with consequent 
ambiguities. A second more phenomenological approach is to try to guess a (string-inspired) effective field 
theory in five-dimensional geometry whose dual (4-dimensional) theory then has QCD-like properties. This 
is known as the AdS / QCD approach by analogy with Maldacena's conjecture, although the term is slightly 
misleading due to the fact that deformed Anti-de-Sitter spaces are usually considered. There are then no 
ambiguities due to choice of compactification manifold, but the choice of 5-dimensional geometry one starts 
with is itself undetermined. Thus, it is important to classify the properties of various geometries and to 
examine the constraints which may be imposed to rule out various alternative geometries. 

In order to determine whether or not a particular geometry (which may or may not have a specific theo- 
retical motivation) approximates QCD well, it is useful to consider more than one geometry and compare 
results coming from the AdS / QCD correspondence in a systematic manner. The aim of this paper is 
to undertake such an investigation, using the static interquark potential as an observable to constrain the 
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theory. The usefulness of this quantity stems from the fact that one can easily separate the regime where 
the semi-classical string calculation is expected to hold (large distances) from the regime where the AdS / 
QCD approach is expected to break down (small distances). 

A five-dimensional holographic hadron model was introduced in [2J. A class of models exhibiting linear 
confinement was examined in [3J H] . In [5] the heavy quark-antiquark potential was considered in a particular 
geometry belonging to this class, and found to be consistent with the Cornell potential [6]: 

V(r) =V + K r-~ + t, (1) 

where the coefficients can be found by fitting e.g. to Charmonium spectra. The original Cornell potential did 
not include the term in r~ 2 , which strictly speaking is not physical. However, such a term is often introduced 
when parameterising lattice determinations of the potential, in order to enhance the fit to data by modelling 
running coupling effects etc. For a review of lattice measurements of the potential, see [7] . Absence of stringy 
corrections in the AdS / QCD calculation corresponds to a lack of quark dynamics. Thus in comparing the 
calculated potential with lattice results, one should use quenched data in which the quarks are not dynamical. 

It has been shown previously that Cornell-like behaviour arises from the AdS / CFT correspondence sub- 
ject to quite general conditions on the geometry [8]. In this paper we consider various geometries in the 
five-dimensional framework, including the recently calculated geometry of [S], in which the back-reaction of 
the 5-dimensional quark condensate fields are used to deform an initial AdS space. It is hoped that such 
a geometry better models the physics on the field theory side of the correspondence, which motivates a 
comparison with other possible geometries. We use a x 2 goodness of fit criterion to compare the potentials 
with lattice data. It is found that no geometry is able to describe the whole range of lattice data for the 
potential (which is to be expected), but that certain geometries are better able to describe the moderate 
and large distance data than the back-reacted geometry of [9]. However, this latter geometry seems to en- 
compass more physics than a similar geometry with the same number of variable parameters. A geometry 
based on a simple quadratic warp factor [5] appears by inspection to agree most closely with the data, but 
this conclusion is argued to be somewhat misleading. 

The paper is laid out as follows. In section [2J we recall how to calculate the heavy quark-antiquark potential 
from a string world-sheet in 5-dimensions spanning a Wilson line on the 4-dimensional boundary [101 [5] . 
In section 3 we present explicit results for the potential using the back-reacted geometry of [9], using the 
parameters fixed as in that paper and the procedure of [5] . A discrepancy is found between the potentials at 
small r values, and in section 4 we consider a more general and systematic procedure for constraining AdS 
/ QCD geometries. Section 5 contains a discussion of the results, and conclusions. 

2 The Heavy Quark- Antiquark Potential 

The short-distance r _1 behaviour of the interquark potential follows from the asymptotic AdS behaviour of 
the bulk space [T], and linear confinement at large distances was demonstrated on quite general grounds in 
[8]. For completeness and to set up our notation, we recall here how to calculate the interquark potential. 

In field theory, the potential between two sources of the gauge field is calculated from the expectation value 
of a Wilson loop W(C) connecting them in space-time. In fact, one has: 

< W(C) >oc e" Ty M (2) 

where the limit T — > oo is understood and V(r) is the separation-dependent potential. According to the AdS 
/ CFT correspondence, the expectation value of a Wilson loop in the boundary conformal field theory is 
related to the extremal area of a string world sheet in the higher-dimensional theory, which spans the curve 
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C [TD]. One has: 

< W(C) >oc e~ A , (3) 

where A is the area of the string worldsheet. Combining this with equation @, the potential (up to a 
constant) is given by: 

V(r) = -1— f d 2 ^detg nm d a X™d p X"\ (4) 
Ana 1 J 

where the Nambu-Goto action for the string worldsheet has been used, and it is understood that this 
worldsheet has extremal area. Here a' is the string tension, g nm the space-time metric, X n (£) are the space- 
time coordinates and £ = (£o, £i) are the worldsheet coordinates. One can now find the form of the potential 
by assuming a given background metric g n m{X), and it is therefore at this point that the model-dependence 
enters. This method was applied in the framework of the AdS / QCD approach in OUT]. Here we assume 
a general form of the metric: 

dz 2 

ds 2 = f(z)dx"dx fl + — , (5) 

where z is the fifth dimension and x^ the 4-dimensional coordinates. If f(z) = z~ 2 , this gives a Euclidean AdS 
space. The possible form of the deformation function f(z) will be considered later. Substituting equation 
(0) into equation ((4|) and choosing £o = t and £i = x (possible because of reparameterisation invariance of 
the worldsheet), one obtains following the method of [S]: 



f z ° 1 
r(zo) = I dx = 2 dz — = 

Jo V/ 



and: 



Vr(zo) = 



z ° dz 
z 2 



f 1- 



4P 



(6) 



(7) 



where z is the value of z at x — 0, r(zo) the interquark separation and g = 1/a' . We have also introduced 
f(z) — f(z)/z 2 , which aids numerical convergence of the integrals. The potential given by equation ([7]) 
diverges and must be regularised [10], as seen by the fact that f(z) — > 1 as z — > from the asymptotic 
AdS behaviour of the metric. Replacing the lower limit of the integral in equation [7] by z — 5, one finds a 
regularised energy: 



V R (zo) = lim 

o — »-U 
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Under general conditions for f(z) the potential will be Cornell-like [ISllH], with asymptotic limits: 
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r — > oo, 



(9) 
(10) 



where p = r 2 (l/4)/(27r) 3 / 2 and {z*, /*} are the values of zq and f(zo) at which the interquark separation 
diverges to infinity. This will occur according to the criterion: 
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which gives the equation for Zq~. 
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= 2fS, 



(11) 
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In general one is only able to solve this equation numerically for Zq , although we will see in the next section 
that an analytical solution is possible for the back-reacted geometry of [9]. It is furthermore useful to 
introduce the following ansatz for the deformation function f(z): 

f(z)^e 2A ^. (13) 



3 Results in a Back Reacted Geometry 

The holographic hadron model of [2] features a five dimensional field dual to the 4-dimensional bilinear quark 
operator qq permeating a fixed 5-dimensional (AdS) geometry. This is an approximation which ignores the 
fact that the action for the 5-dimensional field theory in a curved space-time couples the quark field to the 
metric. Thus, the presence of a bulk field will itself deform the AdS geometry, leading to the deformation 
factor f(z) discussed above. This effect was taken into account in [9], where a differential equation is given 
for the "warp factor" A(z) of equation (fT3|) in terms of the bulk scalar field @. One solves this equation 
subject to a suitable ansatz for this field, and more than one example is given in [9]. We use as our example 
the following back-reacted warp- factor (denoted "Model 1" in [9]): 

2 2 

A(z) = - log z+^z 2 -^z 4 + —z 6 , (14) 
y 1 B 48 32 48 ' V ' 

coming from the following ansatz [2] for the bulk field X dual to the bilinear quark operatoJl: 

X = ^Lz+ a -z\ (15) 

Care is needed to obtain equation (fl4|) from the result given in [9] due to the use of a Euclidean signature 
here. This technicality is explained in appendix |A"1 Using equations (fl"2l [T^|) . we can find the value zq — Zq 
at which the interquark separation diverges. It is given by the equation: 

^(z*) 6 -^(z*) 4 + ^(z *) 2 -l = 0, (16) 
which is a cubic equation for (zq) 2 , whose real solution for general m q and a is: 



216cr 4 — m^er 3 



3a 2 



(17) 



In [5] the values of m q and a were constrained by considering the number of quark colours, meson masses 
and the pion decay constant, thus leading to the values m q = 2.4MeV, a — (326MeV) 3 . Then one has: 

zl = 4.341 GeV~\ / * = 1.395 (18) 



for the extremal values of the parameter z and the deformation function. From equation (|10p . one has for 
these values of m q and a: 

V R -> (0.01178 5 )r, r -> oo (19) 

where g has yet to be fixed. In 5 , this parameter is fixed from the Cornell potential. The lattice parame- 
terisation of [12] gives k = 0.183 in equation ((TJ. Adopting this procedure here, one finds g = 15.55. This 
is much larger than the value g = 0.94 given in [5] (albeit for a slightly different geometry), this difference 
mostly arising from the parameters of the warp factor in the geometry. If one fixes the values of m q and a 
from [9], one has no choice but to fix g in this way in order to obtain the correct long-distance behaviour. 

2 Note that the dilaton (f> is assumed constant in this analysis, consistent with the semi-classical approximation. 

3 This is not to be confused with the string space-time coordinates introduced earlier. 
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Figure 1: The potential Vr(j-) arising from the geometry (|14[) with parameters as discussed in the text 
(solid), alongside the quenched lattice data of [TH [T31 [H] . The constant term has been adjusted. Also shown 
is the Cornell potential obtained from fitting the parameters (dashed), as described in section |U 



The potential can be directly compared with lattice data in the quenched approximation, where dynamics 
of the quarks is suppressed. We use the dat£@ of [HI US H6], shown in figure 3 as a function of r/ro, where 
ro = 0.5fm = 2.5GeV, and with a constant added such that V(ro) = 0. The constant term of the AdS poten- 
tial is also adjusted so that equality is reached with the lattice potential at asymptotically large distances, in 
order to more easily facilitate a comparison. One sees a marked deviation between the two potentials away 
from the asymptotically long-distance regime. Indeed the value of g obtained gives, via equation ([9]), the 
short-distance behaviour: 

Vh(r)--— , r^O (20) 
r 

in marked contrast to the lattice potential used to parameterise the data in [12 , which has e = 0.321 
in equation ([T]). This difference has arisen from fixing m q and a according to [5], then using the long- 
distance behaviour of the Cornell potential to fix the only remaining free parameter g. There are then no 
free parameters left which influence the short-distance behaviour - g is the only parameter governing this 
regime. Although one does not expect the semi-classical theory used here to give a correct description at 
short distances, the possibility exists of using the parameters in the deformation function f(z) to tune the 
long-distance behaviour of the potential. Then g can be used so as to extend agreement between the AdS 
/ QCD and Cornell potentials into the moderate r region, and the parameters m q and a varied to tune the 
long-distance behaviour. This exercise can be repeated for different candidate geometries, and the results 
used to constrain geometrical parameters, and also compare the phenomenological validity of the various 
models. This is the subject of the following section. 

4 Constraining AdS / QCD Models 

Figure Q] compares the potential obtained by fixing the parameters according to [9] with lattice data. This 
gives no quantitative indication of how good the agreement is. To measure the goodness of fit of the 

4 We are grateful to Gunnar Bali for providing files of the data from |12l 113] . Further precise lattice data exists at low r 
values 1 141 115). but this will not affect the results presented in this paper. 
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Figure 2: Variation of the minimum \ 2 /A^.o.f. with the data cut r cut for (a) the back-reacted geometry of 
[H] ; (b) the quadratic warp factor geometry of [5] ; (c) the geometry denoted III in table H] (d) the geometry 
denoted IV in table 0] 



theoretical potential to the data, one may define the standard x 2 variable: 

Wdata , T/ AdS/QCD T/latM 2 

X 2 = E f^^-* (2D 

where the sum is over the A^ata data points, and <7, is the absolute error associated with the i th point. We 
also introduce the number of degrees of freedom: 

■ATd.o.f = A^data — ^Vparamj (22) 

where Afparam is the number of free parameters associated with the geometry. When minimised, a % 2 /A^.o.f. 
of approximately 1 denotes a good fit, and here a problem occurs in that it is not possible to achieve a good fit 
to all data points using the potential of equations (JSJ [8]). This is to be expected given the discrepancy at low r 
values noted in the previous section after fixing parameters at high r. One of course expects the AdS / QCD 
calculation to break down away from asymptotically large r values, due to a multitude of reasons. Firstly, 
there are stringy corrections to the semi-classical approximation employed here. Secondly, the type of five- 
dimensional description used in this paper must in principle arise from compactification of a 10-dimensional 
string theory, none of which physics is well-known. The very small r behaviour is accessible in the field the- 
ory description, and one expects modifications of the naive r _1 behaviour due to the running of the coupling. 

With this level of ignorance, the best one can do is to impose a cut r cut on the lattice data, excluding all 
values satisfying r < r cut . Physically, this corresponds to tuning the AdS / QCD geometry parameters 
using data only in the regime where the calculation is expected to hold. For each value of r cut , there is a 
corresponding minimum value of x 2 that can be achieved by varying the geometry. One then gains some idea 
of where the string calculation breaks down by looking at how the minimum goodness of fit per degree of 
freedom, x 2 ( r cut)/A / d.o.f., varies with r cut . This is shown in figure[2]for the back-reacted geometry considered 
in the previous section, marked as (a) in the figure. One sees that one can fit the theoretical potential to 
the lattice data well for r values down to approximately 4.4 GeV -1 , before the goodness of fit per degree of 
freedom rises sharply as more data is added. As well as the parameters of the warp factor, the coupling g 
is also allowed to vary, together with a constant term added to the potential as this has no effect on dynamics. 



G 



The potential obtained by fitting for a value of r cut such that the minimum \ 2 /-^d.o.f. — 1 is shown in figure 
[TJ One sees that the fit has not improved the description at small distances - indeed, the fitted potential 
actually lies slightly further from the lattice data at low r. Instead the slope at large r has varied in order to 
increase the goodness of fit. This is a clear indication that the string calculation breaks down at moderate 
r values. Of course, there is no guarantee that this geometry corresponds to QCD. In the above analysis 
one is, after all, fitting a function known to be asymptotically linear with some data which also displays 
such behaviour. It is then only natural that at some point the theory and data must diverge, and it is not 
clear that this is really due to the breakdown of the AdS / QCD approximation. However, further insight is 
gained by repeating the above analysis with different geometries. 

We use for comparison the metric of [5] which has a warp factor containing a quadratic term in the expo- 
nent; also a geometry with a two-parameter exponent involving terms in z 2 and z 4 ; finally a three parameter 
exponent involving terms up to and including z e . The different metrics are summarised in table 1, with 
the latter two denoted by III and IV. The first thing to note from figure [2] is that if one takes the criterion 



Geometry 


A(z) 


Parameters 


Back Reacted [9] 
Andreev et. al. [5] 
III 
IV 


m t r 2 ™ q a 4 a 2 6 
48 ■* 32 ^ "r" 48 A 

-z 2 

k lZ 2 \k 2 Z 4 

kiz 2 + k 2 z 4 + k 3 z e 


m q , cr, g, const. 

c, g, const. 
ki, fe, g, const. 
ki, k-2, ka, g, const. 



Table 1: The geometries used to compare AdS / QCD predictions with the lattice data, where A(z) = 
A(z) + log(z) denotes the polynomial terms in the exponent of the warp factor. The coupling g is also 
allowed to vary, as is an arbitrary constant term added to the potential. 

X 2 /Vd. .f. ^ 1 very seriously, it is not possible to distinguish which geometry if any best describes the lattice 
data. Each satisfies this goodness of fit only for r cu t > 4.4GeV _1 . However it is clear that whilst the x 2 
increases sharply for the back-reacted and III geometries as r cut decreases below this value, it is still possible 
to fit the data well using the quadratic warp factor geometry of [5]. This is despite the fact that the back- 
reacted geometry has an additional parameter. Nevertheless, the back-reacted geometry does outperform 
geometry III, which has the same number of parameters but not the same dependence on z in the exponent 
of the warp factor. Thus, it seems that the back-reacted geometry is indeed to some extent modelling the 
physics more correctly than a naive expansion in z containing the same number of parameters. 

An improved fit is obtained by adding an additional parameter, as in geometry IV. Here the warp factor has 
the same polynomial as the back-reacted geometry, but the coefficients are allowed to vary independently. 
From figure [2j one sees that this improves the range over which the lattice data can be fit and that this 
range roughly agrees with the result of the quadratic warp factor. 

For each geometry, one obtains from figure [5] an approximate value r cut = r* ut at which the x 2 P er degree 
of freedom is approximately unity. The 'best fit' parameter values are then the values arising from fits with 
these cuts. These are collected in table HI One sees that when the coefficient of z 4 is allowed to vary inde- 
pendently from the z 6 coefficient, there is no preference for negative values. This may explain some of the 
improvement in the fit with respect to the back-reacted geometry, in which negativity of the z 4 coefficient is 
fixed as a consequence of the Wick rotation from the Euclidean space of [5] together with the interpretation 
of m q and a in terms of quark mass and vacuum condensate terms respectively. 

The potentials arising from each of these 'best fits' are shown in figure El alongside the lattice data. All of 
the potentials produce the correct large r behaviour as expected from the fits and the results of figure [2] 
However, a striking feature of figure [3] is the variation in the behaviour at small r. Geometry III is furthest 
from the lattice data at moderate and small r, which is consistent with the results of figure [2l However, 



7 



Geometry 


* 

'cut 


Parameters 


rScst r it 


Back Reacted 


A QK 

4. .30 


m q 


n (\A Q A 

U.0484 






a 


a A c c n 

0.0550 






const. 


n 9/1/1 
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1 A O 
14. Z 


Andreev et. al. 


A A 

4.U 


c 


A ooo 

0.888 






const. 


A n^;c 
-U.ybo 






g 


l.ZO 


Ill 


6.1 


k-i 


000598 








0.00500 






const. 


0.240 






9 


15.7 


IV 


4.4 




0.00185 






k 2 


0.000500 






fc 3 


0.00100 






const. 


0.252 






9 


5.46 



Table 2: 'Best fit' values for the parameters of table |4] obtained by cutting the lattice data in each case such 
that the minimum X /-^d.o.f. — 1> and then fitting. All parameter units are in terms of GeV. 
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Figure 3: The static interquark potential arising from the 'best fit' examples of each geometry: back-reacted 
geometry (solid); quadratic warp factor (long dashes); III (short dashes) and IV (dot dashed). The lattice 
data is shown with error bars. 
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the potential from the back-reacted geometry does not fare much better. The geometry that appears to lie 
closest to all the data is the quadratic warp factor geometry, appearing to be even better than geometry 
IV which has two more parameters. However, this latter comparison is not meaningful as it is not valid to 
compare the potentials at values of r at which the \ 2 P er degree of freedom is much greater than unity. Using 
f cut = 4GeV from the quadratic geometry (table |3J , this means one is not able to meaningfully compare the 
results for r/r$ < 1.6 in figure [3) Thus even though by eye the quadratic geometry seems to fit the data over 
the whole range of r very closely in figure [3j appearances are in this case deceptive. The fit is not good over 
this range when quantitatively measured using the \ 2 > an d given that the quadratic geometry and geometry 
IV share a similar range over which they agree with the data at high r one should regard the separation of 
the corresponding potential curves in figure [3] as being some measure of the theoretical uncertainty in the 
calculation. 

5 Discussion 

In this paper we have considered the Cornell-like potential arising from the 'Model I' geometry obtained by 
[9] in their analysis of back-reaction effects. Cornell behaviour naturally arises consistent with previously 
found conditions [8], and one can solve exactly for the asymptotic limits of the potential in terms of the 
parameters of this geometry. We examine the potential for the choice of parameters given in that paper. 
Once these are fixed, one can fix the only remaining parameter g by matching to the long-distance behaviour 
of the Cornell potential. This gives a marked deviation between the AdS and Cornell potentials away from 
this regime, with a coefficient of the r _1 behaviour at short distances that shows little agreement with the 
Cornell result. This leads to the more general question of how one should constrain AdS / QCD geometries, 
and we then give an example of a systematic procedure to test the validity of proposed models. 

Two types of question arise. Firstly: of a given set of candidate geometries, which best approximates QCD? 
Secondly: if one knows that a particular type of geometry does indeed represent QCD, where does the AdS 
/ QCD calculation break down? We have investigated both these questions by allowing the parameters of a 
number of geometries to vary simultaneously, and then quantifying the goodness of fit to lattice data using 
a standard \ 2 minimisation. 

The issue of where the string calculation is expected to break down can be investigated by varying the cut on 
the lattice data that is applied to filter out the region contaminated by unknown effects. A given geometry 
should not be trusted in the region where the lattice data increases the \ 2 P er data point to greater than 
approximately unity. For each of the potentials considered, this amounts to r values less than a few GeV -1 . 

Comparing the variation of \ /A^.o.f. with r cut for different geometries, as in figure [21 gives some informa- 
tion on which geometries are closer to QCD. For example, the fact that the back-reacted geometry is better 
able to describe the lattice data (in the regime where we trust the calculation) than a geometry with the 
same number of parameters, is a clear indication that the physics is better described. This is encouraging, 
given that the back-reacted geometry has a better theoretical motivation. However, the simple quadratic 
geometry of [5] gives better results. The fact that the resulting potential is closer to the lattice results at 
small r (see figure [3j may nevertheless be misleading given that the x 2 per degree of freedom is no good in 
this regime, and there is a large theoretical uncertainty associated with the prediction. One can get closer to 
the data from the potential predicted by the back-reacted geometry by allowing more free parameters, with 
the same caveats at low r. Indeed, the fact that adding higher order terms in the warp factor brings the 
potential closer to the lattice data suggests that the good agreement from the quadratic warp factor may be 
a coincidence. 

Another potential problem with the results presented here arises from the fact that we have only included the 
intcrquark potential in the fits, and it is not strictly correct to vary all the parameters of the geometry with- 
out also checking consistency with other observables that can be calculated from the theory. In principle, one 
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Figure 4: Variation of \ 2 /N^.o.i. with the data cut r cut for (a) geometry IV with unconstrained parameters; 
(b) geometry IV with parameter hi constrained to be negative; (c) linear fit to the potential. 



should extend the \ 2 definition of equation pi"]) to include meson masses, decay constants etc. that also de- 
pend on the free parameters. However, the fact that the best fit parameters from the back-reacted geometry 
do not stray too far from the values originally fixed in [9] means that one can be reasonably confident the ad- 
ditional constraints can be satisfied. This is not necessarily true for geometry IV, which has only been fit here 
to the lattice potential. Other observables are not as clean as the lattice potential in that it can be difficult 
to separate the regime where the AdS / QCD calculation is expected to hold from that where it breaks down. 

It may also be the case that varying the parameters over all possible values is unphysical. For example, 
in the back reacted geometry positivity of m q and a follows from their interpretation in terms of a heavy 
quark mass and vacuum condensate of the bilinear quark operator. This then leads to a negative z term in 
the warp factor of the back-reacted geometry. No such constraint is present for the geometry IV, and it is 
therefore interesting to repeat the minimisation of section [4] whilst ensuring that the coefficient A; 2 of the z A 
term in the warp factor remains negative. The resulting variation of x 2 /-^d.o.f. with r cut is shown in figure 
IH One sees that when ki is required to be negative, the resulting x 2 /iVd.o.f. does diverge for a slightly larger 
'"cut value but this is not a significant effect. 

Also shown in figure 2] is the goodness of fit of a straight line to the lattice data, which then gives some idea 
of where the linear regime breaks down. Comparing figures [2] and [4] suggests that the back-reacted geometry 
begins to fail describing the lattice data even for r values still in the linear regime. This is not the case for 
geometry IV and suggests that, whilst the back-reacted geometry seems to include more physics than the 
case of geometry III, the description is still incomplete. 

Of course, the back-reacted geometry adopted here is not unique but relies upon a particular ansatz for the 
quark operator. A higher order expansion for this operator would result in higher order terms in the warp 
factor. It would be very interesting to see if these higher order terms improved the results shown her^|. One 
could also add extra terms to the warp factor which are unconstrained by back-reaction effects and compare 
the resulting fit with the back-reacted case. The resulting phenomenologically determined potential would 
then act as a constraint on various proposed theoretical models. 



5 Recently, the back-reacted geometry of [9] has been extended to include the effect of the strange quark |17| . This does not 
affect the results of this paper, as the nature of the polynomial in the warp factor is unchanged. 
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To conclude, we have used a y 2 goodness of fit measure to compare predictions for the static interquark 
potential from various geometries in the AdS / QCD correspondence with lattice data. It is possible to 
tell the difference between the predictions, with some clearly more favoured than others. A quadratic warp 
factor is favoured by the data, as is a potential having a higher order polynomial warp factor with parameters 
unconstrained by back-reaction. However, it still may be the case that higher order back-reacted geometries 
prove more successful in approximating QCD. 
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A The Back-Reacted Geometry in Euclidean Space 

In [9] back reaction effects are considered starting from the geometry: 

ds 2 = e- 2A( - y Uxadx" - dy 2 , 



(23) 



where we denote the warp factor with a tilde to avoid confusion with A(z) as defined in this paper. Intro- 
ducing z — e y , one then has: 

ds 2 = e- 2A WdXudx* - (24) 
z 1 

Note that instead of using a different symbol for the warp factor in terms of z, we clarify any ambiguity by 
denoting the arguments of this function explicitly To obtain a Euclidean signature one can Wick rotate: 



to obtain: 



ds 2 



dx^ — > idx^, z 



-2AM (^2 + dx 2 



IZ 



d^_ 

r 2 



The "Model I" warp factor of [9] is given in Minkowski space by: 



(25) 
(26) 



Irs i r s m Q 7 m Q <J 4 a fi 

A(z) = log(z) H q -z 2 H 2_z 4 -| z *> 

w BW 48 32 48 



Substituting this in equation (|26[) gives: 



|ds 2 



1 m q 2 m 9°" 4 

— exp — -z — z 

z 2 y \ 24 16 



(27) 



(28) 



We note that the requirement of a positive sign in a warp factor consisting only of A(z) ~ z 2 in Euclidean 
space was noted already in ^4] based on the need for a discrete spectrum of meson masses. 
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